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Indistinguishable particles in two dimensions can be characterized by anyonic quantum statistics more 
general than those of bosons or fermions. Such anyons emerge as quasiparticles in fractional quantum 
Hall states and certain frustrated quantum magnets. Quantum liquids of anyons exhibit degenerate 
ground states where the degeneracy depends on the topology of the underlying surface. Here we present 
a novel type of continuous quantum phase transition in such anyonic quantum liquids that is driven 
by quantum fluctuations of topology. The critical state connecting two anyonic liquids on surfaces with 
different topologies is reminiscent of the notion of a 'quantum foam' with fluctuations on all length scales. 
This exotic quantum phase transition arises in a microscopic model of interacting anyons for which we 
present an exact solution in a linear geometry. We introduce an intuitive physical picture of this model 
that unifies string nets and loop gases, and provide a simple description of topological quantum phases 
and their phase transitions. 

Phases of matter can exhibit a vast variety of ordered states that typically arise from spontaneous symmetry breaking and can 
be described by a local order parameter A more elusive form of order known as 'topological order' ' reveals itself through 
the appearance of robust ground-state degeneracies, but cannot be described in terms of a local order parameter Examples of 
such topological quantum liquids are the fractional quantum Hall states^ where the ground-state degeneracy depends on the 
number of 'antidots' which can be viewed as punctures (holes) in the two-dimensional surface populated by the quant urn Hall 
liquid^. It has long been proposed that topological quantum liquids also occur in certain frustrated quantum magnetd^EEI^I^, 
but it has only been in recent years that strong candidate materials have emerged'" ". While quantum Hall liquids break time- 
reversal symmetry, the exotic ground states of frustrated quantum magnets are expected to preserve time-reversal symmetry. 
As a consequence of this symmetry many unexplored phenomena may appear, including the intriguing possibiUty of topology 
driven quantum phase transitions which is the central aspect of this manuscript. 

In this manuscript we develop an intuitive physical picture for the emerging low energy physics of topological quantum liquids 
and their phase transitions in terms of surfaces and their topology. We thereby provide a visualization of the underlying quantum 
physics, which is in one-to-one correspondence to a detailed analytical framework. Here we consider systems that preserve 
time-reversal symmetry which in this picture will be described by quantum liquids on closed surfaces. Such liquids exhibit 
ground-state degeneracies that depend (exponentially) on the genus of the surface. A section of an extended high-genus surface 
formed by a triangular arrangements of 'holes' is shown in Fig. [T] Through every such hole there can be a flux of the liquid 
populating the surface. An exponential degeneracy then arises from the possible flux assignments through the holes. While in 
the presence of a flux a hole cannot be contracted, we can eliminate the hole in the absence of flux without changing the state of 
the topological liquid. If there is no flux through any of the holes, they can all be removed, and the state of the quantum liquid is 
identical to that on two separated sheets, as shown on the left side in Fig.[T] It is this state that exhibits topological order. On the 
other hand, if there is no flux through the tubes in the interior of the surface (centered around the black lines in Fig. |2]i, we can 
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FIG. 1 : Phase transition in two dimensions. Two-dimensional surfaces with different topologies that are populated by anyonic quantum liquids. 
A quantum phase transition driven by fluctuations of the surface topology connects the anyonic liquid on two separated sheets (on the left) and 
decoupled spheres (on the right). 
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pinch them off. The resuhing state of the quantum Hquid is then identical to that of disconnected spheres, as shown on the right 
side in Fig.[T] This state has neither ground-state degeneracy nor topological order. 

Here we will introduce a microscopic model which energetically favors the absence of flux through the holes or tubes, thus 
dynamically implementing the two topology changing processes mentioned above. The competition of the two processes drives 
a quantum phase transition between the two extreme states. Our model is defined on the 'skeleton' that surrounds the holes in the 
interior of the surface as illustrated in Fig.|2]where the skeleton forms a honeycomb lattice. The fluxes in the tubes are associated 
with discrete degrees of freedom on the edges of the skeleton lattice, corresponding to anyonic particles of the quantum liquicH^. 
The set of degenerate ground states of the liquid is now in one-to-one correspondence with all labelings of the edges consistent 
with a given set of constraints, characteristic to the underlying quantum liquid. ^^^^ 

As a simple example, we consider a quantum liquid of so-called Fibonacci anyoni^lMllS] Hqj-q there are only two possible 
labelings, namely the trivial particle 1 and the Fibonacci anyon r. At any trivalent vertex of the skeleton lattice, there is a 
constraint forbidding the appearance of only a single r-anyon on the three edges connected to the vertex, allowing the following 
possibilities: 

Due to this constraint the edges occupied by a r-anyon form a closed, trivalent net known as a 'string net^^l One might as 
well identify the two degrees of freedom (1, r) with the two states of a spin- 1/2 (1, .[) and thus the same states can be viewed 
as representing the ground states of a Hamiltonian with three-spin interactions enforcing the vertex constraint above (no single 
j-spin around a vertexj^. 

Returning to our model, we can now specify its microscopic terms 

edges e plaqucttcs p 

The first term favors a trivial label £{e) = 1 on the edge e corresponding to the no-flux state. The second term favors the no-flux 
state <j){p) — 1 for the plaquette p. When expressed in terms of the labels £{e), the plaquette flux (j>{p) is a complicated, but 
local expression involving the twelve edges connected to the vertices surrounding a plaquette, see Fig. [2] and is explicitly given 
in the supplementary material. In the absence of the first term ( Jg = 0) the plaquette term will effectively close all holes, and 
the ground state of the above Hamiltonian describes that of the quantum liquid on two parallel sheets as illustrated in Fig. [T] 
The latter is precisely the string-net model first introduced by Levin and Wen^, which is also closely related to another model of 
string nets discussed recently by Fendley Similarly, in the absence of the plaquette term, Jp — 0, the edge term with coupling 
constant ,/e will close off all the 'tubes' thus leading to the ground state of the quantum liquid on multiple disconnected spheres 
as illustrated in Fig. [T] This edge term acts as a string tension in the string net model, or as a magnetic field in its spin model 
representation. 

In the presence of both terms in the Hamiltonian, quantum fluctuations are introduced which correspond to fluctuations of 
the surface. These fluctuations are virtual processes where plaquettes or tubes close off and open depending on the flux through 
them. We can visualize these fluctuations as local changes to the genus of the surface. If the two terms in the Hamiltonian 
become comparable in strength, the competition between the two drives a quantum phase transition between the two extremal 
topologies (see Fig. 1). At this quantum phase transition the fluctuations of the surface become critical and the topology of 
the surface fluctuates on all length scales. We can visualize the (imaginary) time-evolution of this quantum critical state as a 





FIG. 2: Microscopic model. Our microscopic model energetically favors the flux-free states for the holes and tubes (shaded) of the illustrated 
two-dimensional surfaces. For the surface with a triangular arrangement of holes shown in the left panel the anyonic degrees of freedom in 
our model are associated with the edges of the honeycomb lattice skeleton that surrounds the holes in the interior of the surface. For the linear 
geometry of holes on the right the skeleton lattice forms a ladder geometry. 
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FIG. 3: Phase transition in one dimension. Illustration of the quantum phase transition driven by fluctuations of the surface topology in a 
linear geometry which connects the extreme limits of a 'single cylinder' (on the left) and 'two cylinders' (on the right). 

'foam' in space-time, which is reminiscent of the noti on of a quantum foam introduced by John Wheeler for fluctuations of 3+1 
dimensional Minkowski space at the Planck scale'** 

To understand the nature of this transition, we first focus on the linear geometry shown in Fig. |2] In this geometry the 
Hamiltonian becomes 

H = -Jr ^ '5£(r)^i - Jp ^ ^<l>{p)A 1 (2) 
rungs r plaquettes p 

where the first term now only acts on the rungs between the holes (i.e. on those edges of the skeleton which separate two 
neighboring plaquettes), in analogy to the original model. This model exhibits a continuous quantum phase transition between 
the two extreme topologies shown in Fig. |3] This continuous transition is driven by fluctuations of topology. It turns out 
that the gapless theory describing this transition can be solved exactly as discussed in more detail below and explicitly in the 
supplementary material. 

The two extreme topologies connected by this transition in the linear geometry are as follows: In the limit of a vanishing rung 
term, Jr ~ 0, the ground state is that of an anyonic quantum liquid on a single cylinder where all the plaquettes are closed, as 
shown on the left in Fig. [3] For Fibonacci anyons this ground state is two-fold degenerate, with either a r-flux or no flux through 
the cylinder. In the opposite limit of vanishing plaquette term, Jp = 0, we can close off all the rungs and the ladder splits into 
two separate cylinders with a four-fold ground state degeneracy (either a r-flux or no flux in either of the cylinders), as shown 
on the right in Fig. |3] 

In both limits excitations above these ground states are gapped quasiparticles with a gap of Jp or J^, respectively. The first 
excited state above the 'single cylinder' ground state is a r-flux threading a single plaquette, which prevents it from being closed 
as illustrated in Fig.|4^). In the opposite limit of the 'two cylinder' ground state the first excited state is a r-flux through one of 
the rungs, leaving this rung connecting the two cylinders as shown in Fig.j^). Turning on a small coupling 7^ 0, or Jp 7^ 
respectively, these excitations delocalize, but remain gapped and form bands in the energy spectrum, as explicitly displayed 
in Fig. |5^). For large couplings, some of these excitations proliferate and their gap vanishes at the quantum phase transition 
mentioned above. 

The full phase diagram is shown in Fig.|5j)), where we parameterize the two couplings on a circle as Jp = cos 6 and J^ — sin 9. 
Positive (negative) coupling constants indicate that the no-flux (r-flux) states are energetically favored and the two extreme limits 
discussed above then correspond to the points 9 = Q and 9 = 7r/2 on the circle. The continuous phase transition between these 
two distinct topologies occurs for equal positive coupling strengths Jr = Jp, which corresponds to the point 9 = 7r/4 on the 
circle. 

We can visualize this critical point as a quasi one-dimensional quantum foam, with topology fluctuations of the surface on all 
length scales. As a first step, we have performed a detailed numerical analysis of this critical point using exact diagonalization 
of systems with up to 36 anyons. The continuous nature of the phase transition reveals itself in a linear energy-momentum 
dispersion relation, which is indicative of conformal invariance. A detailed analysis of the energy spectrum further allows to 




FIG. 4: Excitations. Plaquette (left) and rung (right) excitations above the two extreme ground states illustrated in Fig. [5] 
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FIG. 5: Left: Energy spectra. Energy spectra of our microscopic model near the decoupling point (6 = 7r/2). The rung excitations shown in 
Fig.|4]form a gapped quasiparticle band well below a continuum of states (shaded). Open symbols show results from exact diagonalization 
of systems with 24 to 36 anyons. These bands are well described by second order perturbation theory around the decoupling point shown as 
solid lines. Right: Phase diagram. The phase diagram of our microscopic model (j2j where the couplings are parametrized as Jp = cos 9 and 
Jr = smO. The gapped topological phases are indicated by the shaded regions. The topology driven quantum phase transition occurs at the 
exactly solvable critical point 9 = 7r/4. An extended critical phase is found in the region 9 £ (tt, 37r/2) around the second solvable (critical) 
point 9 = 57r/4. 



uniquely identify the corresponding conformal field theory (CFT), which in this case turns out to be the 7th member- of the 
famous series of so-called unitary minimal CFTs^^ with central charge c 14/15. This particular identification of a conformal 
field theory is part of a broader scheme which connects the gapless theory of the topology driven phase transition with the nature 
of the underlying anyonic liquid. In the present case of a quantum liquid of Fibonacci anyons we can make an explicit connection 
between the (total) quantum dimension of the anyonic liquid and the central charge of the conformal field theory. 

In fact, the Hamiltonian at this point is even exactly solvable. The key insight leading to this exact, analytical solution is the 
observation that the Hamiltonian of our topological model can be mapped precisely onto a particular version of the restricted- 
solid-on-solid (RSOS) model, which is exactly integrable and directly leads to the above-mentioned CFT^^. This mapping 
expUcitly connects the Hamiltonian at this critical point with an integrable Hamiltonian defined by the Dynkin diagram Dg 




Here the particular labeling of the Dynkin diagram arises from the underlying topological structure of our model. Specifically 
the labels describe the topological fluxes in the two extreme limits of the model as illustrated in Fig. [3] with the limit of a 'single 
cylinder' in the picture on the left corresponding to the blue circles in the Dynkin diagram and the limit of the 'two cylinders' 
pictured on the right corresponding to the green circles. This underlying structure also gives risJ^Sto a representation of the 
Temperley-Lieb algebra^^ which is characterized by the total quantum dimension d — \J2 + of the anyonic liquid, where 
= (1 + v^)/2 is the golden ratio. A more detailed discussion of the exact solution is given in the methods section and the 
supplementary material. 

Varying the couplings in our Hamiltonian there is another way of connecting the two phases depicted in Fig. |3] which is to 
change the sign of both couplings in the Hamiltonian. For opposite sign the two terms now favor r-fluxes through rungs and 
plaquettes, respectively, which again leads to a competition. Interestingly, we find that this competition results in an extended, 
critical phase separating the two topologically distinct phases, as depicted in the phase diagram of Fig.|5] For the full extent of 
this critical phase we again have topology fluctuations on all length scales. However, the gapless theory describing this phase 
turns out to be in a different universality class as compared to the critical point discussed above. These results can again be 
obtained through a combination of numerical and exact analytical arguments, which are detailed in the supplementary material. 
In particular, there is another integrable point in this extended critical phase for equal coupling strengths = Jp, which 
corresponds to the angle Q = 57r/4 in the phase diagram of Fig. |5] and is thus located exactly opposite of the one discussed 
above. Following a similar route one can map the Hamiltonian at this second integrable point to another variant of the RSOS 
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model associated with the Dynkin diagram Dq. The gapless theory at this point then turns out to be exactly the parafermion 
CFT with central charge c — 7/5. The stability of this gapless theory away from the integrable point is due to an additional 
symmetry of our modei^EI Numerically, we find that it extends throughout the whole region where both couplings favor the 
T-flux states all the way to the points 6 = tt and 9 = 3tt/2, where there is no longer a competition of the two terms of the 
Hamiltonian and the ground states have fluxes either through all plaquettes or rungs, respectively. 

Returning to our original discussion of the model (1) on the surface in Fig.[T| the question arises whether we can understand 
the nature of the quantum phase transition here as well. We can explicitly address this question in the context of another kind of 
anyons, the so-called semions-**. Again, there are two possible labelings, the trivial particle 1 and the semion s. The constraint 
now only allows zero or two semion particles s at any trivalent vertex. The set of edges carrying a semion s form loops instead 
of nets and give rise to what is known as a 'loop gas'*^^. In its spin- 1/2 representation (where f, | now stand for 1 and s) this 
model is known as the honeycomb version of the toric cod^, where the string tension corresponds to a magnetic field. This 
model exhibits a continuous quantum phase transition in the 3D Ising universality class""^^' with topology fluctuations on all 
length scales. Mapping the 2+1 dimensional semion system to its three-dimensional classical counterpart, the quantum foam 
then corresponds to the critical fluctuations of domain walls in a 3D Ising model at its critical point. For other kinds of anyonic 
liquids, the nature of the topology changing transition is in general unknown and remains an intriguing open problem with the 
possibility of new universality classes. For a liquid of Fibonacci anyons there has been a recent discussion of quantum critical 
behavior b y F endley from the perspective of ground-state wavefunctions and their respective correlators in terms of conformal 
field theor>P3. 

Finally, in order to explore the broader context of our models we complete our analysis by considering the complete set of 
possible excitations present in these models. An excitation different from the ones already discussed arises when relaxing the 
constraint which for every trivalent vertex of the skeleton lattice forbids the occurrence of a single r-flux. If we allow for this 
possibility, we are left with a r-flux entering the vertex through one tube, but not leaving it through another tube in the skeleton 
plane as illustrated in Fig.|2] Instead we can think of the remaining r-flux at such a vertex as leaving through one of the liquid 
sheets surrounding the skeleton lattice. This piercing of the liquid by a r-flux corresponds to a vortex excitation of the liquid 
and is illustrated as a 'chimney' in Fig. |6] These vortex excitations break time-reversal symmetry and turn out to all possess the 
same chirality (indicated by the red arrow in Fig. [6]l. This is only possible if the anyonic liquid on a given sheet itself possess 
a given chirality. Since the entire system exhibits time-reversal symmetry, this means that the two anyonic liquids on the two 
sheets must have opposite chirality. Vortices associated with chimneys on opposite sheets thus also have opposite chirality as 
illustrated in Fig.|6] (In fact, a vortex in one sheet can be related to a vortex in the opposite sheet by dragging a vortex through 
a 'hole' connecting the two sheets. Moreover, we can create a 'hole' connecting the sheets by glueing together two vortex 
excitations on opposite sheets.) This conceptual perspective of two anyonic liquids with opposite chirality giving rise to a time- 
reversal invariant model connects with and allows a visualization of a more abstract mathematical description of these models, 
namely doubled non-Abelian Chern-Simons theorieJ^H. It remains an intriguing question to formulate our topology driven phase 
transitions within such a field theoretical framework. 

In this manuscript we have developed a general, unifying framework to formulate topological aspects of quantum states of 
matter for systems preserving time-reversal symmetry and of their phase transitions. In a simple and intuitive picture they 
are described in terms of fluctuations of two-dimensional surfaces and their topology changes. Our framework gives a new 
perspective on how to broadly discuss quantum phase transitions out of topologically ordered states of matter, which has so far 
been largely unexplored territory due to the lack of a local order parameter description amenable to a Landau-Ginzburg- Wilson 
theory. This description of time-reversal invariant anyonic quantum liquids is also expected to advance our understanding of 




FIG. 6: Excitations of the anyonic liquid. Vortex excitations of the liquid indicated by the 'chimneys' posses a chiral edge mode. 
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spin liquid states and their phase transitions in recently proposed materials of frustrated quantum magnetism and other strongly 
correlated systems. Our unifying perspective on string nets and quantum loop gases will also allow to construct a large variety 
of new microscopic models for topological phases. 

Acknowledgments - We thank M. Freedman, X.-G. Wen, and P. Fendley for stimulating discussions. Our numerical simula- 
tions were based on the ALPS libraries^. A. W. W. L. was supported, in part, by NSF DMR-0706140. 



Methods 



Identification of conformal field tlieories. To characterize the conformal field theory (CFT) of the critical points in the linear 
(ladder) geometry, we rescale and match the finite-size energy spectra obtained numerically by exact diagonalization for systems 
with up to L = 36 anyons to the form of the spectrum of a CFT, 



E — El L - 
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+ h+h 



(4) 



where the velocity v is an overall scale factor, and c is the central charge of the CFT. The scaling dimensions h+h take the form 
h = + n,h ~ + n, with n and n non-negative integers, and and hP are the holomorphic and antiholomorphic conformal 
weights of primary fields in a given CFT with central charge c. The momenta (in units 2n/L) are such that = h — h or 
kx = h — h + L/2. Using this procedure we find that for the critical point at 6' = 7r/4 the rescaled energy spectrum matches the 
assignments (28 1 of the 7th member'^" of the famous series of so-called unitary minimal CFTs^^ with central charge c = 14/15. 
Similarly, at the point 9 — 57r/4 we find the rescaled energy spectrum to match that of the parafermion CFT with central 
charge c— 7/5. For the calculated energy spectra and the details of these assignments we refer to the supplementary material. 



Exact analytical solution. The Hamiltonian in Eq. (|2i can be solved exactly for interaction strengths corresponding to angles 
9 — it/ A and 9 = 57r/4 in the phase diagram of Fig. 5] This exact, analytical solution of the gapless theories at these points 
unambiguously demonstrates the continuous nature of the related quantum phase transitions and points to generalizations of 
these gapless theories for other kinds of anyonic liquids. The key observation underlying this exact solution is the emergence 
of the Dg Dynkin diagram from the topology of the surface associated with the ladder model as depicted in the right part of 
Fig.|2] Each labeling of the edges of the 'skeleton' graph which corresponds to that surface, denotes one of the states spanning 
the Hilbert space of the system. A crucial step is to consider a different 'pants decomposition'^^ of this surface and to perform 
a basis change to a new basis which corresponds to the labeling of the skeleton lattice of this alternative pants decomposition. 
Explicitly, this basis transformation can be written as 
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Here, (-Fb"d)a' denotes the so-called F-matrix, which is a generalization of the familiar 6j symbols of angular momentum 
coupling in conventional quantum mechanics and is known for any anyonic liquid^^. Note that associated with the even- 
numbered indices of these labels, which correspond to the original rung labels c, on the right, there is the flux through the 
cross-section of the surface on the left, denoted by a label di = \ or di — r. Similarly, associated with the odd-numbered 
indices, which correspond to the original plaquettes, there is a pair of fluxes through the two cross-sections of the surface at 
the position of the plaquette on the left, denoted by a pair of labels, {ui, hi). This pair of labels can assume four values, i.e., 
(fli, bi) = {(1, 1); (r, 1); (1, r); (r, r)}. The (fusion) constraints at the vertices where the labels (a.;, fe;) and di±i meet then turn 
out to be precisely the condition that they be adjacent nodes on the Dq Dynkin diagram of Eq. (|3]l. For example, a local label 
(fli, bi) = (r, t) at an odd-numbered index i only allows for labels di-i = 1 and di^i = t at the neighboring even-numbered 
indices. This is reflected in the Dynkin diagram by the appearance of a line that connects the label (r, r) to both labels 1 and r. 
The importance of the just described basis change consists in the fact that in the new basis the rung and plaquette terms, i/[ and 
Hf, respectively, of our ladder Hamiltonian 



H = -Jr ^ i/[ - Jp Yl , 

i even i odd 



(5) 
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where 




D HI 
DHf 



for i even, 
for i odd. 



(7) 



The characteristic 'D-isotopy' parameter of this Temperley-Lieb algebra, D = \J\ + Lp^ = 2 cos (tt/ 10), is precisely the total 
quantum dimension of the underlying Fibonacci anyon liquid. We have thereby established a remarkable, explicit connection 
of the one parameter of this emerging algebraic structure, the 'D-isotopy' parameter of this Temperley-Lieb algebra, and the 
single most characteristic parameter of the underlying anyonic liquid, namely its total quantum dimension. This observation 
points to a generalization of such a connection for other quantum liquids. Written in this form, the resulting Hamiltonian for the 
Fibonacci anyon liquid turns out to be precisely that of the (integrable) restricted solid-on-solid (RSOS) statistical mechanics 
lattice model based on the Dg-Dynkin diag ranP, as obtained in the standard fashion from the transfer matrix of the RSOS 
lattice model. For further details we refer to the supplementary material. 

General framework. We have explicitly formulated the concept of a topology driven phase quantum phase transition mainly 
in the context of a single anyon theory, namely that of Fibonacci anyons. However, we note that these concepts apply in great 
generality to any anyon theory in which there can be an arbitrary number of anyons subject to a set of fusion rules / constraints. 
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Supplementary Material 



I. FIBONACCI ANYONS 



The d egree s of freedom in our microscopic models are so-called Fibonacci anyons, one of the simplest types of non-abelian 
anyons'^ The Fibonacci theory has two distinct particles, the trivial state 1 and the Fibonacci anyon r, which can be thought 
of as a generalization (or more precisely a 'truncated version') of an 'angular momentum' when viewing the Fibonacci theory 
as a certain deformatiorP^ of SU(2). We will now make this notion more precise and illustrate it in detail. In analogy to the 
ordinary angular momentum coupling rules, we can write down a set of 'fusion rules' for the anyonic degrees of freedom which 
are analogs of the Clebsch-Gordon rules for coupling of ordinary angular momenta, 

1X1 = 1 IXT — T — TXl T X T — I + T , (8) 



where the last fusion rule reveals what is known as the non-abelian character of the Fibonacci anyon: Two Fibonacci anyons t 
can fuse to either the trivial particle or to another Fibonacci anyon. In more mathematical terms, these fusion rules can also be 
expressed by means of so-called fusion matrices Nj whose entries {Nj)j]^ equal to one if and only if the fusion of anyons of 
types ji and j2 into j is possible. The fusion rules are related to the so-called 'quantum dimensions' dj of the anyonic particles 
by 



N,\d,) = d,\d,) , 



(9) 



where \dj) is the ('Perron Frobenius') eigenvector corresponding to the largest positive eigenvalue of the 2 x 2-matrix Nj. [The 



sense in which these numbers are 'dimensions' will become apparent in section 11 A 1 below.] For the particles in the Fibonacci 
theory the quantum dimensions are di = 1 and dr ^ (p = {I + a/5)/2 and the total quantum dimension of the theory is then 
given by = (^^. dp 1/2 = yTT^. 

To define our Hamiltonian, some additional indegredients of the theory of anyons are required. In analogy to the 6j'-symbols 
for ordinary SU(2) spins, there exists a basis transformation F that relates the two differents ways three anyons can fuse to a 
fourth anyon, depicted as 





(10) 



The left hand side (l.h.s.) represents the quantum state that arises when anyon a first fuses with anyon b into an anyon of type 
e, which, subsequently, fuses with anyon c into an anyon of type d. Similarly, the right hand side (r.h.s.) denotes the quantum 
states that arises when anyon b first fuses with anyon c into anyon type / which, in turn, fuses with anyon a into anyon type d. 
Whilst keeping all external labels, the types of the three anyons (a, b, c) as well as the resulting anyon type d fixed, the states on 



the l.h.s and r.h.s. are fully specified by the labels e and /, respectively. Eq. ( 10 1 says the so-specified states are linearly related 
to each other by the so-called F-matrix-'^ with matrix elements {F^^^^l. 

In general, the F-matrix is uniquely defined (up to 'gauge transformations') by the fusion rules through a consistency relation 
called the pentagon equatiorP^ Similarly, the braiding properties of anyons are given by the so-called i?-matrix (which however 
is not needed here) that is uniquely determined by the hexagon equation^^. 

For the Fibonacci theory, it is straightforward to verify that in most cases there is only one term on the right-hand-side in 



Eq. ( 10 1, e.g. by choosing two or three out of the four anyons a, b, c, d to be r-anyons. For these cases the consistency with the 
pentagon and hexagon relations then yields that the corresponding F-matrix elements equal to 1. There is only one configuration 
that gives rise to F-matrix elements that are non-trivial: If all anyons are r-anyons, e.g. a — b — c— d — t, both the 1- and the 
T-fusion channels appear, and the F-matrix takes the explicit form 

As a final ingredient to explicitly derive our Hamiltonian, we have to introduce the so-called modular S'-matrix that relates the 
anyon "flux" of species b through an anyon loop of species a to the case without anyon loop by 



Ob 

5J 



(12) 
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For the case of Fibonacci anyons, the S'-matrix takes the expUcit form 



SI S; J D\'f -1 

There is an important relationship between the modular S'-matrix and the matrix encoding the fusion rules, introduced in the 
paragraph above Eq. (j9]l: the modular S'-matrix diagonalizes the fusion rules, the 'Verlinde Formula', 

St, {Nat, S^c - Si (14) 



(repeated indices are summed) where S^ denotes the adjoint of the unitary matrix S. The eigenvalues of the matrix (Na) are 
thus -5T, and the largest (positive) eigenvalue, the quantum dimension da, can be seen to be 

da=^- (15) 

Due to the unitarity of the modular S'-matrix one immediately checks that the total quantum dimension equals 

D=ji- (16) 



II. THE LADDER MODEL 



In this section we will discuss details of the "ladder model" in a one-dimensional geometry, whose Hamiltonian is given by 
Eq. (2) in the main part of the paper. We start by defining the Hamiltonian in detail, and then discuss the gapped topological 
phases, critical phases, and the exact solutions. 



A. The Hamiltonian 

1. Explicit expression 

To establish a notation for the basis states we consider the skeleton lattice inside the high-genus ladder geometry as shown in 
Fig. [7] The basis states are given by all admissible labeling of the edges of the skeleton with 1 or r particles, subject to the vertex 
constraints given by the fusion rules. The number of basis states, B^, of the ladder with L plaquettes and periodic boundary 
conditions is given by 

{ai,bi,Ci} {«!,.. .11,} 

where Ni are the fusion matrices of Fibonacci theory as introduced above. The largest eigenvalue of the matrix Ni is the quantum 
dimension di. Thus, the leading behavior of the traces for large L is, 

L 

Bl^ {d^nd^....d,,Y = Y,{d\)'^-\dlf^{\ + ^^f^D^^. (18) 

{ii...i_L} fe=0 

The Hilbert space thus grows asymptotically, for large L, as a power of the square of the total quantum dimension D^. 
The Hamiltonian (as given in Eq. (2) of the main part of the paper) 

H = -Jr ^ ^£(r),l - Jp Y^ ^0(p),l (19) 
rungs r plaq p 

consists of two non-commuting terms, the rung term which is diagonal in the chosen basis, and the plaquette term which depends 
on the four edges of the plaquette p, and the four adjoining edges. By inserting an additional anyon loop of type s into the center 
of the plaquette, we can project onto the flux through this additional loop (and hence the flux through the plaquette) by the 
following procedure (for a derivation see the following subsection) 
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FIG. 7: The high-genus surface in a ladder geometry and the 'skeleton' of the fusion graph that defines the Hilbert space. The trivial particle 
1 or the Fibonacci anyon r can occupy the links of the ladder skeleton, subject to the vertex constraints given by the fusion rules of Fibonacci 
anyons. 
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The additional s-loop is inserted by performing a sequence of F-transformations: 
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(20) 



(21) 
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(22) 



(23) 



The ladder geometry has a local duality between the inside and outside: the inside of the rungs is dual to the plaquettes. The 
only difference is that the rungs connect two different cylinders, while the plaquettes connect the same space (the "outside"). 
The duality can be made exact by using "twisted" boundary conditions where the ends of the ladder are connected according to 
Oi — bL+i and 61 = a^+i (so that the ladder looks like a Moebius strip). Indeed, our exact diagonalization results confirm that 
the excitation spectra are identical under exchange of the couplings and Jp for twisted boundary conditions. However, in the 
case of periodic boundary conditions (ai = a^+i, &i — 6l+i), which we shall focus on in the following, this duality is only up 
to degeneracies. 
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2. Bigger (mathematical) picture 



So far our discussion in this 'Supplementary Material' has been largely focused on detailed algebraic manipulations. In this 
subsection we wish to give a brief idea of the general bigger picture of topological field theories which underlies these detailed 
manipulations. At the same time we will provide a deeper understanding of the so-called 'Levin-Wen model' within this context. 

In the main text we have given a physically motivated description of the Levin-Wen model in Figs. 1 and 2, leading to the 
Hamiltonian in Eq. (1) of the main text. Let us now give a more abstract description of it. The most general Levin- Wen 
Hamiltonian has two kinds of terms: the vertex type (not discussed so-far as a term of the Hamiltonian) and plaquette type. Let 
us consider a surface S, and a trivalent graph F (which we called 'skeleton'in the main text) embedded in that surface. (The 
sole role of the surface E, which in the leftmost picture of Fig. 1 of the main text is just a parallel plane sitting in between the 
two depicted sheets, is to give a well defined meaning to the notion of a 'plaquette'; namely, all complimentary regions of Y\T, 
i.e. the complimentary regions of the graph T within the surface S, are plaquettes.) We always enforce strictly the condition 
that three labels meeting at a vertex must satisfy the fusion rule. (This is another way of saying that we have set the coupling 
constant of the 'vertex term' in the most general Levin- Wen Hamiltonian to infinity.) As a result, we obtain a Hilbert space 
called L{T, S) consisting of the Hilbert space spanned by all admissible labelings of the trivalent graph F: a labeling of F is an 
assignment of a label in a label set Ic to each edge of the graph, [Iq — {1, r} in the previous subsection|l|, and the labeling is 
admissible if the three labels around each vertex satisfy the fusion rules. 

Now, there exists another vector space, which brings about the connection with the actual surfaces that were drawn in Figs. 1 
and 2 of the main text. In particular, when C denotes a so-called modular category (for a precise definition, which we do not 
need at the moment, see e.g. Ref. l4TT l which basically denotes a theory of 'anyons' and their corresponding 'fusion rules' such 
as the one described in the previous subsection |l] then the vector space L{T, S) is the same as a Hilbert space Vc{Sr) (for a 
definition see e.g. Ref. 1411 1 of an associated Topological Quantum Field Theory (TQFT) corresponding to the 'modular category' 
C: specifically let Nr be the thickening of the graph F to a handle-body (drawing a cylinder around each edge of the graph), and 
Sr be the boundary surface of iVr, then L{T, E) = Vc{Sr)- In the language of TQFT, any 'pants-decomposition' of the surface 
Sr is known to lead to a basis of Vc {Sr), which corresponds to the vector space spanned all possible fusions of the labelings on 
F. 

This interpretation of the Hilbert space L{T, E) gives rise to a transparent derivation of the plaquette term, Eq. (23 i, in the 
Levin-Wen model. To derive this expression, we use the identification of L{T, E) with Vc(S'r)- The cth row of the modular 
5-matrix of the modular category C can be used to construct a projector lUc that projects out the particle with a label c through 
a plaquette. In other words total flux c through a plaquette p can be enforced by inserting ujc into a plaquette p. The projector 
turns out to be of the form 



(24) 



where [a] denotes a loop labeled by a as the one drawn in Eq. ( 12 1. In order to see that this performs the task let us insert a flux 



with label b thought the loop [a], resulting in the figure drawn on the l.h.s. of Eq. ( 12 1, which we denote in symbols by [a]{b) 



When we now perform the sum in Eq. ( 24 1 we obtain, upon making use of Eq. ( 12 1 




_ J_ ca 
^ D i—ia "^c 



I 



(25) 
(26) 



15IT6I 



where we have used the unitarity (plus reality and symmetry) of the modular S-matrix, as well as Eq.s 

Therefore, the plaquette term (^^(p) i is implemented by inserting the projector lui = ■ [a] into the plaquette p. Now the 

detailed steps leading to Eq. (23 1 are easy to understand: The insertion of uii into the plaquette is written explicitly in Eq. (20i. 
In the subsequent equation, first an i^-move is applied to the two lines connected by the dotted line, and subsequently four more 
F moves counterclockwise around p are implemented as drawn; finally removing the resulting bubble, we obtain the explicit 
form of the plaquette term written in Eq. ( [23] l. 

The mathematical context for the Levin- Wen model is the Drinfeld center Z{C) or quantum double of a unitary fusion category 
C. The label set Ic for the Levin- Wen Hamiltonian is the isomorphism classes of simple objects of C. It is known that a unitary 
fusion category is always spherical. By a theorem of M. Mugei^'^, the Drinfeld center of any spherical category is always 
modular. It follows that the Drinfeld center of any unitary fusion category is always modular Moreover, if the spherical category 
C itself is modular, then Z{C) is isomorphic to the direct product of the conjugate C* and C, where C* is obtained from C by 
complex conjugating all data. Our main example is one of those special cases, where C is the Fibonacci theory. 



12 



0-- -o- -o- -o- 



FIG. 8: Possible basis configurations in the presence of one plaquette excitation. The two ground states (without a flux through the plaquette) 
and three excited states (with a r-flux through the plaquette) are linear combinations of these basis states. Solid lines denote r-anyons while 
dashed lines symbolize the trivial particle 1. 



The decomposition of Z{C) hints directly at the appearance of Dynkin diagram at the critical point in one-dimensional 
geometry: indeed, the two phases connected by the critical point are the Fibonacci theory and the doubled Fibonacci theory with 
label sets {1, r}, and {(1, 1), (1, r), (r, 1), (r, r)}, respectively. Based on this, it is natural to expect that the two sets of fusion 
rules will fit together in a compatible way at the critical point, which is nicely illustrated by the structure of the Dg Dynkin 



diagram in Fig. 17 which underlies the exact solution of this critical point (Section 11 D below) 



B. Topological phases 

We start the detailed discussion of the phase diagram with the two distinct gapped non-abelian topological phases: the 'single 
torus' phase where all plaquettes are closed at 6* = and the 'two tori' phase with closed off rungs at 9 = tt/2. A finite-size 
scaling analysis of the splitting of the ground state degeneracies and the energy gap shows that the phases extend over a wide 
range of parameter space as illustrated in the phase diagram (Fig. 5b of the main part of the paper). In this section we discuss 
the low-lying excited states in these phases, give their explicit wave functions at the exactly solvable points, and a perturbative 
expansion for their dispersion away from these points. 

1. The 'single torus' phase at — 7r/2 < 9 < n/4 

To describe the lowest excited states we consider the trivially solvable point 9 = where Jr = 0. In the ground state there is 
no flux through any of the plaquettes, and they all can be closed, thus reducing the high-genus ladder to a single torus (see Fig. 
3 in the main part of the paper). There are two degenerate ground-states configurations with either no flux or a r-flux through 
this torus. 

Similarly, we can deduce the degree of degeneracy for the lowest excited state by considering the topology of this state. In 
the lowest excited state, one plaquette flux is present which yields the reduced topology (as compared to the high-genus ladder) 
and the associated skeleton shown in Fig. 4b of the main part of the paper. Closing all but one plaquette this skeleton allows for 
5 different 1 , t coverings, illustrated schematically in Fig. |8] In order to obtain the anyon-fluxes through the excited plaquette, 
a basis transformation (consisting of a F- and a S'-transformation) of the reduced basis is performed which yields that there are 
three r-fluxes through each plaquette. Thus, the lowest excited state at 6* = is 3i-fold degenerate. Tuning away from 9 = 
these 3L excitations delocaUze and form a three-fold degenerate band. 

2. The 'two tori ' phase a/7r/4<5<7r 

At the point 9 = tt/2 (trivially solvable) the ground state has no r-anyons on the rungs of the ladder. The rungs can hence be 
cut which yields an effective topology of two separate tori. Of the four degenerate ground states three are symmetric and one is 
antisymmetric under j/-reflection. The lowest excited state is a r-anyon flux through a single rung. The fusion rules then require 
a flux through both of the two tori, and this state is hence only L-fold degenerate. Tuning away from 9 = tt/2, these states 
delocalize into a non-degenerate band. 

3. Perturbation expansion for the quasiparticle bands 

Over a broad range of parameters the quasi-particle excitations are well described (see Fig. 5b of the main part of the paper) 
by a second order perturbative expansion around 9 = tt/2, with a dispersion given by 

AE{Jp, Jr, fcx) = Jr-^ cos(A:,) - -^[1 + 2cos(fc,)] - ^^|^2cos(2fc,). (27) 
Due to duality, this result equally applies for coupling parameters 9 close to 9 = 0, with Jr and Jp interchanged. 
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C. Gapless theories 

In this section, we discuss the critical points (9 = 7r/4, 57r/4) and the extended critical phase in the ladder model. We first 
discuss the gapless theories in terms of numerical results and then present analytical arguments leading to an exact solution for 
the two critical points (9 = tt/4, 5tt/4). 



1. Critical point at = tt / 4 
( numerical findings from exact diagonalization) 

At equal positive values of the two coupling constants (Jp = Jr,9 — tt/4:), the system has a linear energy-momentum 
disperson relation with the finite-size spacing between energy levels vanishing linearly in 1/L. This indicates that the two 
adjacent, gapped topological phases (Fig. 5) are separated by a continuous quantum phase transition and a critical point that is 
described by a 2D conformal field theory (CFT). To characterize this CFT, we rescale and match the finite-size energy spectra 
obtained numerically by exact diagonalization for systems with up to i = 36 anyons to the form of the spectrum of a CFT, 

E ^ EiL+ — (-— + h+h] , (28) 



L V 12 

where the velocity v is an overall scale factor, and c is the central charge of the CFT. The scaling dimensions h+h take the form 
h = + n, h = Jv^ + n, with n and n non-negative integers, and and hP are the holomorphic and antiholomorphic conformal 
weights of primary fields in a given CFT with central charge c. The momenta (in units 27r/L) are such that = h — h or 
kx ~ h — h + L/2. Using this procedure, we find that for the critical point at 6' = 7r/4 the rescaled energy spectrum matches the 



assignments (28 1 of part of the Kac-Table of the m ~ 9 unitary Virasoro minimal CFT of central charge c = 14/15, as shown 
in Fig. [9] In Fig. 10 we list all relevant primary fields of this CFT which appear and their corresponding scaling dimensions. 



It turns out that only the Kac-Table primary fields (pr.s with s = odd appear, and those with s = 5 have multiplicity two (the 
associated states on the ladder being in the bonding/antibonding sectors of 'transverse momenta' ky = 0, vr), all others having 
multiplicity one. These are precisely those Kac -table primary fields which occur in the so-called (D,A)-modular invariants] of 
the m = 9th Virasoro minimal CFT of central charge c~ 14/15. 

To illustrate how the ground-state degeneracy changes at this critical point from a two-fold degeneracy for the 'single cylinder' 
limit (Jr = 0) to a four-fold degeneracy for the 'two cylinders '-limit (Jp = 0), we can follow the evolution of eigenenergies 
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FIG. 9: Exact diagonalization: Energy spectrum at the critical point (9 — 
it/4) for a ladder with L — 12 holes and 36 anyons. The energies have been 
rescaled so that the two lowest eigenvalues match the CFT scaling dimensions. 
The open boxes indicate the primary fields of the 7th minimal model with 
central charge c = 14/15. The topological symmetry sectors are indicated 
with symbols 1 = yi,i, t = yr,T and r + 1 = yi,T. 



FIG. 10: CFT fields: Scaling dimensions h(^r,s) + /i{r,s) 
of the primary fields in the {D, A) modular invariant of 
the 7th minimal model with central charge c = 14/15. 
On the right, we give momentum and topological sym- 
metry assignments of these primary fields for our micro- 
scopic model. 
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FIG. 11: The energies of the lowest lying energy states around the critical point {Q — 7r/4) as a function of the 'dimerization' 6. Results are 
shown for system size L = 10. 



in the vicinity of this critical point as shown in Fig. 11 Moving away from the critical point (9 = 7r/4) corresponds to a 
dimerization of the model; in an alternative basis choice, discussed in detail in section 11 D it becomes apparent that the rung 
and plaquette terms alternatingly act on even and odd 'sites'. For 6 \ tt/A, the four-fold ground-state degeneracy is lifted 
with one of the four ground states approaching the field with rescaled energy 2/45 (ky — 0), and two degenerate ground states 
moving to a rescaled energy 2/15 {ky = and ky — tt). The single first excited state in this gapped phase softens towards the 
rescaled energy 4/15 at the critical point. As we move into the adjacent gapped phase for 9 < tt/A only the field with rescaled 
eigenenergy 2/45 moves back towards the ground-state, while the two other fields move upwards in energy and form a three-fold 
degenerate excited state. 



2. Extended critical phase for 6 G (tt, Stt, 2) 
(numerical findings from exact diagonalization) 

For negative coupling parameters Jp,Jr < 0, we find an extended critical phase around the point of equal coupling strength 
which in our circle phase diagram is opposite to the critical point discussed above. For the whole extent of this critical phase 
we can match the finite-size energy spectra to the parafermion CFT with central charge c ~ 7/5. This theory is part of 
the sequence of Zj. -parafermion CFTs with conformal weights A^^ — ''''i^o'' ~ T^' where j — 0, ^ 1, fc/2, |m| < j (and 
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and Table 11. 



j — m = integer), in the notation ol'*^. The details of the assignments for A: = 8 can be found in Fig. 

In order to verify that the critical phase around the exactly soluble point 9 = 5tt/4 extends to the vicinity of the decoupling 
points 9 = TT and 9 = 37r/2, we consider an effective model where we fix all rung occupations to r-anyons. This assumption is 
exactly true at the decoupling point 6 = 3tt/2. Implementing this constraint significantly reduces the size of the Hilbert space 
and allows us to numerically study this effective model for larger system sizes with up to 48 anyons. 
The effective Hamiltonian in the reduced Hilbert space is given by 

iJ-ff = -LJr -JpY. ^Hp),i ■ (29) 
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FIG. 12: Exact diagonalization: Energy spectrum at the critical point {9 — 
57r/4) for a ladder with L — 8 holes and 24 anyons. The energies have been 
rescaled so that the two lowest eigenvalues match the GET scaling dimensions. 
The open boxes indicate the primary fields of the Zg-parafermion model with 
central charge c = 7/5. The topological symmetry sectors are indicated with 
symbols 1 = t = yT,T and r + 1 = j/i,^. 



EIG. 13: GET fields: Scaling dimensions + /i(j „,) 
of the primary fields in the Zg-parafermion GET with 
central charge c = 7/5. On the right, we give momentum 
and topological symmetry assignments of these primary 
fields for our microscopic model. 



The first term is a constant, and can thus be omitted which then turns the actual value of Jp irrelevant. A positive Jp corresponds 
to the limit 9 \ 37r/2, while a negative Jp allows to study the limit 9 y 37t/2. 

For positive Jp, we find that the splitting of the ground state degeneracies goes to zero for 1/L 0, and the energy gap 
approaches a finite value as shown in Fig. [14] This further supports the stability of the gapped topological phases up to, but 
excluding, the points 9 ~ n and 9 — in/ 2 in our phase diagram. 

For negative Jp, the rescaled energy spectrum of this effective model is critical and again matches (with much higher accuracy 



than at = tt) the Zg parafermion conformal field theory with central charge c — 7/5 as shown in Fig. 12 We can hence 
conclude that the whole quadrant 9 g (tt, 3tt/2) is occupied by an extended critical phase described by the same conformal field 
theory as the exactly solvable point 9 — 5tt/4. 

Approaching the endpoints of this extended critical phase at 6* = tt and 9 = 3tt/2, the low-energy spectrum collapses into a flat 
band resulting in an extensive ground state degeneracy below an energy gap of size 1 at the points 9 = tt and 9 = 37r/2. Moving 
beyond these 'decoupling points' where one of the terms in the Hamiltonian vanishes, this extensive ground-state degeneracy 
is split again and a gap opens for 6' < tt and 9 > in/ 2, respectively, as the system enters the two gapped, topological phases 
discussed above. 
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FIG. 14: Energy gap AiJ(l/ L) / Jp between the first excited state and the ground state, as well as the splitting of the ground state degeneracy, 
SE{1/L) I Jp, for the effective model Eq.|29] The two ground states become precisely degenerate only in the thermodynamic limit. The results 
indicate that the energy gap extrapolates to a finite value. Since the effective model is valid in the limit 6 3n/2+, the gapped topological 
phase extends all the way up to this point. 
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Topological stability of the critical phases 



Both critical theories have a large number of rescaled energies (28 1 that are smaller than two. These eigenenergies are 
associated with operators whose correlation functions decay with scaling exponent h+h<2. Such operators are relevant in the 
renormalization group sense. This means that any operators O with scaling dimensions (=rescaled energies) h +h < 2 which 
is invariant under all symmetries of the Hamiltonian may appear as an additional term in the latter and can thus drive the system 
out of the critical phase into a gapped phase or a different critical phase. For a critical phase to be stable there must hence exist 
a symmetry in the model such that the identity field (associated with the ground state) belongs to a different symmetry sector 
than all other fields (j> with h+h<2 and — (fields atkxy^O do not obey the translational symmetry of the Hamiltonian). 
Indeed, our model has an additional topological symmetry^^ that can stabilize the critical phases: There can be either no flux 
(denoted as 1-flux) or a r-flux entering the periodic ladder from one side, and a 1- or a r-flux leaving the ladder as illustrated in 



Fig. 16 There are hence three possibilities for possible flux assignments: (i) no flux is entering from above, and no flux is leaving 



[Fig. 16 1], (ii) a r-flux is entering and leaving [Fig. 16 1], or, (iii) a r-flux is entering from one side, and leaves through one or 



several plaquettes as shown in Fig. 16;). For each operator, one of the three scenarios applies and we can explicitly determine 



the topological sectors by considering the following hermitian symmetry operator (which commutes with the Hamiltonian) 

Y\a,b,c) ^ n(^^^)a''^(^^;)^'"1«'^^'.c>, (30) 



where \a, b, c) — \ai, &i, ci, 02, 62, C2, ol, 6l, cl) are labels according to Fig. tI This operator inserts additional r-loops 
parallel to the two 'spines' of the ladder As in the case of the plaquette term Eq. (23 1, this is done by connecting them to the 
ladder with 1-particles. The flux through each of these two additional r-loops can be either 1 or r, where a 1-flux yields a factor 
of SI/ S\ — (p, and a r-flux gives S^/S^ = — f^^ (note that a S'-transformation has to be performed in order to obtain the flux 
through the additional r-loops). Hence there are three possible eigenvalues ofY: yi.i = ip'^ (scenario i), yr.r = f^'^ (scenario 
ii) or T — —ip^^ip — —1 (scenario iii). 



We numerically evaluate the topological symmetry sectors in the two critical phases (see Tables 1 and 11, and Figs. |9 12 



and 15 1. At the critical point separating the topological phases {9 ~ 7r/4), we find that the relevant operators can be classified 
according to s = 1 Vi.i, s = 3,7 <-*■ yr.T, s = 5 ^ Vi.t- In particular, only one operator, (t>{2,i), is in the same 
topological symmetry sector as the ground state, i.e. the identity field 0(1, i). It is this field that drives the system out of the 
critical phase when varying the coupling constant 0. With the scaling dimension of this operator being a; = 2/3 the gap opens 
as AE{9) (X \0 — 7r/4|'' on either side of the critical point, where v = 1/(2 — 2/3) = 3/4. In the second critical phase, 
9 E (tt, 37r/2), the topological symmetry assignments of the relevant operators are given by ?■ = ^ r — 2,6 ^ yT,T, 
r = 4 <-+ ?/i In particular, there is no relevant field in the same topological symmetry sector as the ground state, which implies 
that there is no symmetry-allowed relevant operator in this gapless theory and the critical point must be part of an extended 
gapless phase. This observation demonstrates that our observation (from exact diagonalization studies) that the extended critical 
phase in the quadrant 9 E (tt, 37r/2) is described by the same Zg parafermion CFT with central charge c = 7/5 is correct. 



a) no flux (1-flux) b) r-flux c) (r + l)-flux 




FIG. 16: Topological symmetry sectors: a) No r-flux enters or leaves the ladder, b) A r-flux enters from one side and leaves at the other side, 
c) A r-flux enters from one side and leaves through a plaquette. 
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D. Analytical solution 



Our ladder model defined by the Hamiltonian in Eq. (2) in the main part of the paper can be solved exactly at the two critical 
points 9 = Ti/A and = 57r/4 (see the phase diagram in Fig. 5 of the main text). The key observation leading to this exact 
solution is that the topological structure of our model implies that its Hilbert space is in fact built on the so-called Dg-Dynkin 



diagram, which is drawn below in Fig. 17 



O — 

(1,1) 



-o- 

I 



-o- 




FIG. 17: Dynkin diagram Dg 



The Dynkin diagram indeed appears very naturally: let us make a change of basis for our Hamiltonian as illustrated in Fig. 18 
This new choice of basis (drawn on the left), which arises from a different decomposition of the high-genus surface, is related 
to the original one (drawn on the right) by a simple F-transformation. In particular, consider the new basis in the left part of 



Fig. 18 with the even-numbered 'sites' (which correspond to the original rungs) we associate a label di — 1 ox di — t (the 



flux through that cross-section of the surface). With the odd-numbered 'sites' (which correspond to the original plaquettes) we 
associate a variable consisting of a pair of labels, (oj, hi) which can assume four values, i.e., (oj, hi) = (1, 1), {ai, hi) — (r, 1), 
(fli, hi) = (1, r) and (a^, hi) — (r, r), and denotes the pair of fluxes through the two cross-sections of the surface at the position 
of the plaquette. The allowed fusion channels at the vertices where variables (ci, hi) and di±i meet then correspond precisely 



to the condition that they be adjacent nodes on the Dynkin diagram of the Dq Lie algebra, as illustrated in Fig. 17 above. For 
example, a local label (a^, hi) — (t, t) at an odd-numbered 'site' i allows for labels di-i = 1 and di-i = r at the neighboring 
even-numbered sites, which is reflected in the fact that label (r, r) is connected by a line to both labels 1 and r in the Dynkin 
diagram. 




FIG. 18: Two possible basis choices corresponding to different decompositions of the high-genus surface. The basis drawn on the left is used 
in formulating the exact solution: the rung and the plaquette terms alternatingly act on even or odd 'sites' i. 

In summary, the elements of this new basis of the Hilbert space on which the Hamiltonian acts are of the form 

\d) := |...,ai_i,ai,ai+i,...) , (31) 

where aj [— dj if j is even, and = {aj, bj) if j is odd] denotes a point on the ZJg-Dynkin diagram representing the flux through 
the high-genus surface at the 'site' j of the chain. The sequence of aj must satisfy the condition that a^+i is a nearest neighbor 
site of aj on the Z?g-Dynkin diagram. 

In this new basis, the rung and plaquette terms Hf and Hf of our ladder Hamiltonian 

H=-JpJ2Hf- Jr J2 HP , (32) 

i odd i oven 

take on the following fornJSl 

s=l,T a'. 

^M^)-E(^:a._,a.,,)i.(<^^._,a.,,)d;l<)- (33) 
d'. 



19 



In fact, these terms can be seen to form a representation of the Temperley-Lieb alg which arises from the Dg-Dynkin 
diagram, and has "d-isotopy" parameter D — ^Jl + (p'^ — 2cos(7r/10), the total quantum dimension of our Fibonacci theory. 
Specifically, consider the operators constructed from the components Va — sin(a7r/10) (a = 1, 6, vi — W2 = 

W3 = V{t,t)^ V4 = V(^y vc^ — vg — f(i,r) = ^(t,i)) of the ('Perron Frobenius') eigenvector corresponding to D, the largest 
positive eigenvalue of the adjacency matrix of the ZJg-Dynkin diag 



where 



These operators form a known representatiorP^ of the Temperley-Lieb algebra with "d-isotopy"- parameter D, i.e. 



D e, 



[e„e,]=0 for\i-j\>2. 



(34) 



(35) 



Now one can check that the rung and plaquette terms, Eq. (33 1, of the Hamiltonian in the new basis, Eq. (32 1, are proportional 
to these operators, i.e. 



Hf = ^ foi' * odd' ^ Z) ^* ' even. 



(36) 



The Hamiltonian Eq. (33 1 is in fact that correspon ding to the (integrable) restricted-solid-on-solid (RSOS) statistical me- 
chanics lattice model based on the ZJg-Dynkin diagrarrP^. Specifically, the two-row transfer matrix T :— T2T1 of this lattice 
model 




is written in terms of Boltzmann weights W[i] assigned to a plaquette i of the square lattice 

Ti:=]JW[2n], and T2 := J| W[2n + 1] 



(37) 



with 



sm u 



The parameter it > is a measure of the lattice anisotropy, 1 is the identity operator, and 



(38) 



(39) 



The Hamiltonian of the so-defined lattice model is obtained from its transfer matrix by taking, as usuai^, the extremely 
anisotropic limit, < u <C 1, 



T = exp{-a{B. + ci) + 0{a^)}, a = 



D sin[7r/10] 



«1, 



yielding 



H=-y -ei 
^ D ' 



(40) 
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Since, due to Eq. (|36 1, the operators '-^ ei' are nothing but the rung and plaquette operators, we have thus demonstrated that 
"the RSOS statistical mechanics model based on the Dynkin diagram coincides with the Hamiltonian, 



the Hamiltonian of 
Eq. ( 32 1, of our ladder model. 

The RSOS model based on Dq is knowri^-''^^' to provide an (integrable) lattice realization of the {D, A) modular invarianl^^l 
of the 7th unitary minimal CFT of central charge c = 14/15. In particular, the Hamiltonian of Eq. (2) of the main text at 
angle 9 ^ tt/A will yield the spectrum of that CFT. This exact analytical result is borne out precisely by our numerical (exact 
diagonalization) studies reported in subsection (IIC 1 1. This CFT with central charge c = 14/15 describes the quantum critical 
point of a 1 + 1 D quantum system, our ladder model. While we cannot make an exact statement for the related ^ + 1 D quantum 
model, we note that Fendley has recently discussed a 2 + 1 D quantum critical point from a 2 + D perspectiv^i^by considering 
a one-parameter family of wavefunctions connecting the ground-state wavefunctions of the two extreme limits of the 2 + 1 D 
model ([T]). For a certain value of the parameter he finds a conformal quantum critical point whose ground-state correlators are 
written in terms of this same c = 14/15 CFT. 

Another version of this lattice model yielding in the anisotropic limit the negative, — iJ, of the Hamiltonian in Eq. (40 1 is 
also integrable and provide a lattice realization of the Zg parafermionic CFT of central charge c = 7/5. In particular, the 
Hamiltonian of Eq. (2) of the main text at angle 6* = 57r/4 will yield the spectrum of that CFT. Again, this exact analytical result 
is borne out precisely by our numerical (exact diagonahzation) studies reported in subsection ( IIC2| i. 



III. THE HONEYCOMB LATTICE MODEL 



In this section, we discuss details of the "honeycomb lattice model" whose Hamiltonian is given by Eq. (1) in the paper 
first define the plaquette term of the model and then discuss two limiting phases of the model. 



We 



A. The Hamiltonian 



In analogy to the plaquette term in the ladder model, Eq. (23 i, the plaquette term of the honeycomb lattice model (Eq. (1) in 
the main text) is defined by 



30(p),i 




s=l,r a',/3',-y' 

S',e'X' 



a I b 




where the additional two edges of a plaquette are reflected in two additional F-transfomiations. Again, we can parametrize the 
coupling constants on a circle as Jp = cos(6') and Jg = sin(6'). 



B. Excitations 



We briefly mention the elementary excitations of this model. In the 'two-sheets' phase, which corresponds to couplings = 
( Jp — \, Je — 0), the elementary excitation is a single plaquette with a r-flux giving rise to a single 'hole' as illustrated on 
the left in Fig. [19] These excitations are gapped with a gap size of Jp and will delocalize for small couplings Je 7^ forming 
quasiparticle bands. Similar to the ladder model the dispersion of this quasiparticle band can be calculated perturbatively around 
the 'two-sheets' limit. 

In the opposite limit of 'decoupled spheres', which corresponds to couplings — -k 12 (Jg = 1, Jp = 0), the elementary 
excitation is a 'plaquette ring' where all edges around a given plaquette have r-fluxes, as illustrated on the right in Fig. 19 
Again, a perturbative analysis allows to qualitatively describe the quasiparticle band. 
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